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Abstract —In this paper, joint designs of data routes and 
resource allocations are developed for generic half-duplex multi- 
carrier wireless networks in which each suhcarrier can he reused 
hy multiple links. Two instances are considered. The first instance 
pertains to the general case in which each suhcarrier can he time- 
shared hy multiple links, whereas the second instance pertains 
to a special case in which time-sharing is not allowed and a 
suhcarrier, once assigned to a set of links, is used hy those 
links throughout the signalling interval. Novel frameworks are 
developed to optimize the joint design of data routes, suhcarrier 
schedules and power allocations. These design problems are 
nonconvex and hence difficult to solve. To circumvent this 
difficulty, efficient techniques based on geometric programming 
are developed to obtain locally optimal solutions. Numerical 
results show that the designs developed in both instances yield 
performance that is superior to that of their counterparts in 
which frequency-reuse is not allowed. 

Index Terms —Power control, geometric programming, mono¬ 
mial approximation, time-sharing, self-concordance. 

I. Introduction 

The prospect of having ubiquitous high data-rate wireless 
services is leveraged by the versatility and portability of the 
communication devices that will form the nodes of future 
wireless networks. These devices will be able to perform 
various functions including sending, receiving and/or relaying 
data to other nodes. As such, it is expected that future wireless 
networks will not possess a predetermined topology, but rather 
an ad hoc one that encompasses many existing and upcoming 
network structures including current and relay-aided cellular 
networks m, El. 

Given the stringent limitations on the spectrum available 
for wireless communications, providing high data-rate services 
banks on sharing the spectrum by multiple users which results 
in potentially significant interference. One option to mitigate 
interference is to use the Orthogonal Frequency Division 
Multiple Access (OFDMA) technique, wherein a set of or¬ 
thogonal narrow-band subcarriers are exclusively assigned to 
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each user. This technique offers several advantages including 
design simplicity and resilience to frequency-selective fading. 
In spite of these advantages, rate-effective utilization of the 
available spectrum may require the OFDMA subcarriers to be 
used simultaneously, rather than exclusively, by multiple users. 
This is especially the case when the network is composed of 
essentially separated clusters. In contrast, for tightly coupled 
networks exclusive usage of subcarriers can be more beneficial 
from a rate perspective El, El- 

In order for a wireless network to be able to support the 
reliable communication of high data rates, the scarce resources 
available for the network must be carefully exploited. Such 
resources include the spectrum available for communication, 
time, and the typically low power of the wireless nodes. Proper 
exploitation of these resources involves choosing the optimal 
routes of the data flows, the optimal powers to be allocated 
by the nodes to each subcarrier, and the optimal scheduling 
and possibly duration over which the subcarriers are assigned 
to various links. Although these tasks have traditionally been 
performed separately, they are interrelated and performing 
them in isolation may incur a significant loss in performance. 

To avoid the aforementioned drawback, we devise a joint 
optimization framework that incorporates data routing, sub¬ 
carrier scheduling and power allocation in the design of a 
generic multicarrier network. The network is generic in the 
sense that it possesses an ad hoc topology and its nodes can 
assume multiple roles simultaneously including being sources, 
destinations and/or relays. This framework is centralized, in 
the sense that the design is performed by a central entity that 
is aware of the network parameters. Hence, this framework 
can be seen as a benchmark for distributed and potentially less 
comprehensive designs. In this framework each subcarrier can 
be reused by multiple links, and the nodes acting as relays 
operate in the half-duplex mode, i.e., a node cannot send and 
receive at the same time on the same subcarrier. The objective 
of the design is to maximize a weighted-sum of the rates 
injected and reliably communicated over the network. Weights 
are assumed to be known a priori, but can be adapted over 
time to account for fairness issues and to maintain a desired 
quality of service. 

Two instances of networks are considered. In the first 
instance, each subcarrier can be time-shared by multiple links, 
thereby resulting in continuous subcarrier scheduling vari¬ 
ables. The role of these variables is to determine the fraction of 
time during which a subcarrier is used over a particular link. In 
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contrast with the first instance, in the second one time-sharing 
is not allowed and a subcatTier, once assigned to a set of links, 
is used by those links throughout the signalling interval. This 
instance results in binary subcarrier schedules, which using a 
particular change of variables, are incorporated in the power 
allocation constraints. This results in a design problem that 
is significantly easier than its general counterpart considered 
in the first instance. It is worth noting that the first instance 
is a generalization not only of the second instance, but also 
of instances in which frequency-reuse is not permitted a. 
As such, the framework considered in this instance offers 
significant performance advantages over currently available 
designs, but at the expense of increasing dimensionality and 
design complexity. This instance provides an inherent tradeoff 
between performance and design complexity; the design com¬ 
plexity can be reduced by restricting the number of links that 
can reuse a particular subcarrier. 

The optimization problems arising from the joint design in 
both instances are nonconvex and hence difficult to solve. To 
overcome this difficulty, a logarithmic transformation is used 
to cast the original problem in a form that, for all but a few 
constraints, complies with the geometric programming (GP) 
standard form E). The constraints that are not compatible 
with that form are approximated by monomial expressions 
that correspond to their first order Taylor expansion around 
a given initial point ||71. Using an exponential transformation, 
the resulting approximation can be cast in a convex form. A 
refinement of this approximation can be obtained by iterative 
updating of the initial point. In particular, we use the so- 
called iterative monomial approximation technique, wherein 
the solution of one convex approximation is used as the 
initial point in the following iteration. Under relatively mild 
conditions, this technique is guaranteed to yield a solution 
of the Karush-Kuhn-Tucker (KKT) system corresponding to 
the original problem il. Numerical results show that the 
designs developed in both instances yield performance that 
is considerably superior to that of their counterparts in which 
frequency-reuse is not allowed. 

In comparison with currently available designs, the ones 
presented herein are the first to attempt designing data routes, 
subcarrier schedules and power allocation jointly when the 
subcarriers are both time-shared and frequency-reused. In par¬ 
ticular, the contributions in this paper include; 1) introducing 
the concept of simultaneous time-sharing and frequency-reuse 
of subcarriers; 2) casting the joint design of data routes, 
subcarrier schedules and power allocation in a framework 
that is amenable to GP-based optimization; 3) providing a 
simplified approach that enables mustering a considerable 
portion of the gains offered by the full joint design, but with a 
significantly lower complexity; and 4) developing an efficient 
polynomial complexity algorithm for the special case in which 
the subcarriers can be frequency-reused but not time-shared. 
This work builds on the results obtained in ||9| and ifTOl . How¬ 
ever, the exposition herein is more comprehensive and includes 
additional examples, a simplified approach and complexity 
analysis. 

The paper is organized as follows. Section |II] provides an 
overview of currently available design techniques. Section Hill 


explains the system model and design objective. Section |IV] 
considers the joint design of data routes and power allocations 
when time-sharing of subcarriers is allowed. The comple¬ 
mentary instance in which time-sharing is not allowed is 
addressed in Section lYl The complexity of the proposed 
algorithms are examined in Section |VT] Numerical results are 
provided in Section IVlIl and Section lVIIIl concludes the paper. 
For completeness, the GP standard form and the monomial 
approximation technique are provided in Appendix lAl and the 
derivation of the results pertaining to complexity is provided 
in Appendix Ifil 

II. Related Work 

In this section we provide an overview of the currently 
available techniques for routing and resource allocation in 
wireless networks. A plethora of techniques is available for 
optimizing each aspect in isolation, but significantly fewer 
ones consider their optimization jointly. 

Resource allocation in wireless networks constitutes the 
task of determining the power allocated for each transmis¬ 
sion and the fraction of time over which a particular sub¬ 
carrier is assigned to that transmission. Instances in which 
resource allocation techniques were developed are provided 
in El, m-na for various network scenarios. For instance, 
power allocation techniques for single-carrier cellular systems 
and ad hoc multicarrier systems were developed in El 
and El, respectively. To enable more effective utilization 
of resources, power allocations were optimized jointly with 
binary-constrained subcarrier schedules. For instance, the de¬ 
signs developed in El and El rely on the premise that 
each subcarrier is exclusively used by one node and the 
solutions obtained therein are potentially suboptimal. When 
the binary constraint on the subcarrier schedules is relaxed 
allowing the subcarriers to be time-shared by multiple nodes, 
the optimal power allocations can be shown to be the water¬ 
filling ones Gl; a related problem was considered in El for 
a case in which the nodes experience self-noise. 

Further improvement can be achieved by joint optimization 
of resource allocations and routing Q, El-El- For instance, 
a method for obtaining jointly optimal routes and power 
allocations was developed in El for the case in which 
the nodes were restricted to use orthogonal channels for 
their transmissions. In a complementary fashion, the case in 
which the power allocations are fixed was considered in El- 
Therein, a heuristic was developed for optimizing the data 
routes and subcarrier schedules jointly. 

Capitalizing on the potential gains of incorporating power 
allocation jointly with data routing and subcarrier scheduling, 
the authors considered a generic network in which the nodes 
can assume multiple roles at the same time and each subcarrier 
could be either used exclusively by one link or time-shared 
by multiple links El- Although the designs provided in El 
offer an effective means for exploiting the resources available 
for the network, these designs restrict the subcarriers to be 
used exclusively by only one link at any given time instant. 
Such a restriction may not incur a significant performance 
loss in tightly coupled networks 121, but in networks with 
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clustered structures, this restriction can be quite harmful. For 
unclustered networks, frequency-reuse may result in a substan¬ 
tial increase in the interference levels. However, if properly 
exploited, frequency-reuse can yield valuable performance 
gains. The effect of frequency-reuse was considered in single¬ 
channel networks in m for the case in which the data rates 
are restricted to assume discrete values, and in HD for the 
case in which the nodes use superposition coding. 

In the current work, we will consider the joint optimization 
of power allocations, subcarrier schedules and data routes in 
the design of generic multicarrier networks with frequency- 
reuse, and with and without time-sharing. As such, these 
designs generalize currently available ones, and will subse¬ 
quently offer a significant improvement over their perfor¬ 
mance. 

A summary of this review and a comparison to our work is 
presented in Table J] 


III. Problem Statement 


A. System Model 

We consider a multicarrier wireless network of N nodes, 
each with one transmit and one receive antenna, and a fixed 
power budget, Pn, n € M = {1,2,-•• ,N}. The network 
operates over a frequency-selective broadband channel of 
bandwidth Wq, which is partitioned into K frequency-flat 
narrowband channels, each of bandwidth W — Node 
are assumed to be capable of simultaneously transmitting, 
receiving and relaying data. This assumption is generic, in 
the sense that constraining some nodes to perform a subset 
of tasks can be readily incorporated in the formulations that 
will be developed hereinafter. For tractability, the nodes are 
assumed to always have data ready for transmission M, and 
for practical considerations, the relaying nodes are assumed 
to operate in a multi-hop, rather than cooperative, half-duplex 
mode iflTl . 

The nodes are connected with L wireless links, each com¬ 
posed of K subcarriers and the set of all links is denoted by 
C = {1, 2, • • ■ , L}. The coefficient of the fc-th subcarrier of 
link £ connecting node n to node n' is denoted by the complex 
number which comprises pathloss, shadowing and fading. 
An instance of such a network with N = A nodes and K = 2 
subcarriers is depicted in Figure [T] 

In addition to their desired signals, the nodes receive a su¬ 
perposition of noise and interference due to the transmissions 
of other nodes in the network. Denoting the signals transmitted 
to and received by node n on the fc-th subcarrier by and 
y^n\ respectively, we can write + 

ik') (k) 

, where denotes the set-minus operation and 

denotes the corresponding zero-mean additive Gaussian noise 
with variance Aq. Assuming, as before, that £ S £ is the link 
connecting node n to node n', it can be seen that the signal- 
to-noise-plus-interference ratio (SNIR) observed by node n' 
on subcarrier fc of link I is given by 


SNIR(£, fc) 


Fnfelfc-nri'l 




( 1 ) 


( 2 ) _ ( 2 ) ( 2 ) ( 2 ) 
*2 ~ *3 *4 



Fig. 1: An exemplary network with N = A, K = 2 and D = {1, 2}. 
The objective is to maximize a weighted sum of i.e., the rate 

injected at node n and intended for destination d. 


where pnk is the power allocated by node n to the fc-th subcar¬ 
rier. The second term in the denominator of ([T]i represents the 
aggregate interference observed by node n' on subcarrier fc of 
link £. When the nodes transmit Gaussian distributed signals, 
the maximum data rate that can be reliably communicated on 
this subcarrier is given by IFlog 2 (l -F SNIR(£, fc)). 

For ease of exposition, we divide both the numerator and 
denominator of ([T]i by WNq and we use gtk to denote the 

normalized channel gain, between any two nodes 

n, n' e Af. 


B. Network Topology 

The considered network can be represented by a fully- 
connected weighted directed graph with N vertices and L = 
N{N — 1) links. To facilitate enumeration of links, the link 
from node n to node n' will be labelled by £ = {N — l)(n — 
1) -F n' — 1 if n < n' and by £ = {N — l)(n — 1) -F n' 
if n > n' . The sets of incoming and outgoing links of node 
n G J\f are denoted by C-{n) and £+(n), respectively, and 
the connectivity of this graph can be captured by an incidence 
matrix, A = [ani\, where = 1 if £ G £+(n), Oni = — 1 if 
£ G £-(n) and Oni = 0 otherwise M- 


C. Design Objective 

Let V = {1, ■ ■ ■ ,£?} be the set of all destination nodes, 
where V C A/". Let be the rate of the data stream 
injected into node n G Af and intended for destination 
d G V. The objective of our joint design is to maxi¬ 
mize a weighted-sum of the rates injected into the network, 
i.e., where are non¬ 
negative weights satisfying J2neN\{d} = 

1 . 

Assigning weights to the injected rates provides a conve¬ 
nient means for controlling the quality of service (QoS); a 
higher weight implies a higher priority. Weights are typically 
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TABLE I: Related Work 


Platform 

Routing 

Power Allocation 

Scheduling 

Frequency-reuse 

Reference 

Multiuser system 

X 

/ 

/ 

/ 


Uplink cellular 

X 

/ 

X 

/ 

unn] 

OEDMA 

X 

/ 

/ 

X 

112J, [13], IMl, 115J 

Generic network 

/ 

/ 

X 

X 

IE] 

Mesh network 

/ 

X 

/ 

X 

In] 

Generic network 

/ 

/ 

X 

/ 

II 8 J, [19J 

Generic network 

/ 

/ 

/ 

X 


Generic network 

/ 

/ 

/ 


This work 


assigned a priori, but can be adapted to meet QoS require¬ 
ments M- Varying the weights enables us to determine the 
set of all rates that the proposed design can simultaneously 
achieve. 

Having described the system model, in Section |IV] we 
will characterize the constraints that must be satisfied by the 
routes, the subcarrier schedules, the data rates and the power 
allocations. 

IV. General Case; Routing and Resource 
Allocation with Time-sharing 

We consider the case when each subcarrier can be both 
reused and time-shared by multiple links. This case generalizes 
the cases in which either frequency-reuse or time-sharing of 
subcarriers is not considered, e.g., 0. After characterizing 
the constraints that must be satisfied by the network variables, 
we will formulate the cross-layer design as an optimization 
problem. Unfortunately this problem is nonconvex and to 
obtain a solution of its KKT system, we will use an iterative 
GP-based technique that is guaranteed to converge to such a 
solution. 

A. System Constraints 

In this section, we derive the mathematical constraints that 
must be satisfied by any feasible set of data routes, time¬ 
sharing schedules and power allocations. 

1) Routing Constraints: Let x^^f} be the data flow intended 
for destination d G V on subcarrier k G K, of link I G C. The 
flows, {xg‘1^}, and the injected rates, are related by 

the flow conservation law, which must be satisfied at each 
node. This law stipulates that the sum of flows intended 
for any destination d G V at each node must be equal to 
zero M- Applying this law to the current network and using 
the incidence matrix in Section lTlI-Bl it can be seen that } 
and must satisfy the following constraints: 

anix^f = nGJ\f\{d},dGV. (2) 

£eCkeK 

The flow conservation law implies that the rate of data leaving 
the network at d G V equals the sum of the data rates 
injected into the network and intended for this destination. 
Hence, we can write injected 

rates, are non-negative, and since the network is 

represented by a directed graph, the flows, must be 

also non-negative. Hence, 

s^n'’ > 0> n e A/'\{d}, dGV, (3) 


Ad) 


> 0 , 


G C, k G K, d GV. 


(4) 


2) Scheduling Constraints: Considering both time-sharing 
and frequency-reuse requires introducing a set of variables 
to characterize the fraction of time over which a particular 
subset of links utilize the same subcarrier. To do so, let 
be the fraction of the signalling interval during which 
links ti,... S £ are simultaneously ‘active’ on subcarrier 
k G 1C', the remaining L — m links in £ are ‘silent’ on 
this subcarrier. Without loss of generality, we will write the 
indices in an ascending order, i.e., fi < • ■ • < Im- For 
notational convenience, let T be the set of all the subcarrier 
time-sharing schedules. The cardinality of F is given by 
|r| = (t) = K{2^ — 1). For instance, consider a 

network with £ = 3 links and K = 1 subcarriers. In this case, 
F = 72^3> 71 ^ 2 . 3 } and |r| = 7. To 

see the role of F, consider the schedules in Figure In this 
figure, 7 ^^^^ = 0.5, = 0.2, 7 ^^^ 3 = 0-3^ and all the other 

elements in F are zero. 


Signalling interval fork=l 




7P’(50%) 


71 ^ 3 ( 30 %) 

-> < > 


□ active 

□ silent 


Fig. 2: An exemplary scheduling table for a network with £ = 3, 
K =1. 


Note that the fact that the channels are assumed constant 
over the signalling interval implies that only the time-sharing 
schedules (i.e., entries of F) affect the rate expressions, irre¬ 
spective of the particular time interval over which the subcar¬ 
riers are time-shared. In other words, horizontal displacement 
of the shaded blocks in Figure |2] does not affect the rate 
expressions. 

The number of variables in F grows exponentially with the 
number of links, L. This renders the incorporation of F in the 
joint optimization computationally prohibitive. In most cases 
this complexity can be significantly reduced without incurring 
heavy performance losses. For instance, if the network is 
tightly coupled, high interference levels render the reuse of 
subcarriers on multiple links less beneficial. In such a case, 
restricting the reuse of a subcarrier to a fewer links may 
incur negligible deterioration in performance but reduces the 
number of variables significantly. To take advantage of this 
observation, we limit the number of links that can reuse a 
particular subcarrier to / <C £. By performing this restriction. 
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the number of elements in F is reduced from K (2^ — l) 
to which, for small /, is polynomial in L. For 

instance, if at most two links are allowed to reuse a particular 
subcarrier at any given time, i.e., 1 = 2, the number of 
elements in F reduces to . It is worth noting that lim¬ 

iting the number of simultaneous transmissions, I, inherently 
offers a trade-off between the performance and complexity. In 
particular, as I increases, the available resources are utilized 
more efficiently. However, our simulations suggest that most 
of the gain of time-sharing and frequency-reuse is accrued by 
only considering / < 3 simultaneous transmissions. 

For feasible time-sharing schedules, the elements in F must 
be non-negative and, to ensure no overlapping in time, the 
total time over which any subcarrier fc € /C is used must not 
exceed the length of the signalling interval. These constraints 
imply that 

F > 0, elementwise, (5) 

E E Vfce/C. ( 6 ) 

m=l (i-'-em&C 

Note that summations in (| 6 ]) characterize the number of links, 
TO, that reuse a particular subcarrier k. For instance, for the 
case in which at most two links reuse this subcarrier, the left 
hand side (LHS) of (|6]l can be expressed as 7 ^^^ -|- 

Nodes cannot broadcast data to multiple destinations at the 
same time, that is, at any time instant, node n can have at 
most one active link on subcarrier k. Hence, the time-sharing 
schedules corresponding to multiple outgoing links of node n 
must be zero. This can be represented as 

+ E E = 0 ’ 

”1=3 ^3 '^meC 

£iGA£ 2 e/:\{^i},fcG/C, (7) 

where = max{0, Qni}, that is, = 1 if f G C+{n) and 
zero, otherwise. 

To enforce the half-duplex constraint, we must ensure that 
no two links, ii G C-{n) and £2 G £+(n), can be active 
on the same subcarrier k € K, at the same time. This implies 
that all the time-sharing schedules that correspond to li and 
£ 2 , i-C-, Tij ^ j TO = 2,..., /, must be zero. Since all the 
entries in F are non-negative, these constraints can be written 
as 




i: 


(fe) 


= 0 , 


m—Z 

£2 ^ \ k G /C, ( 8 ) 


where a~^ = | min{0, a„f}|, that is, a~^ = —1 if f G C-{n) 
and zero, otherwise. Note that 0 and ® take effect only 
when ^ 0 and 0 +^ 0 +^, ^ 0 , respectively. 

3) Power Allocation Constraints: To facilitate the design, 
we replace the node power variables {pnk} with link power 


variables {qik\, which are related by the following transfor¬ 
mation; 

Pnk = max qik, n € Af, k G IC. (9) 

iGC+(n) 

To gain a better understanding of the transformation in (|9l), 
we note that 0 implies that, of all the links in £+(n), only 
one element in the set {qik}ieC+(n)^ W n € J\f, k G JC, 
can assume a strictly positive value. Now, 0 indicates that 
this value is the power allocated by node n to subcarrier k. 
Using 0, we will formulate our design in terms of {qgk} 
instead of {pnk}- These variables must satisfy the following 
non-negativity constraints; 

qek >0, ecC, kGfC. (10) 

In a practical network, the nodes are likely to have indi¬ 
vidual power budgets which bounds the total power used by 
each node on all subcarriers. To capture this constraint, we 
note that only the subcarriers scheduled to outgoing links 
contribute to the power consumption of each node. More 
specifically, if G £+(n), then all the time-sharing schedules 
that correspond to ii contribute to the power consumption at 
node n. This constraint can be written as 

keK£i^C+(n) 7n=2e2---tmeC 

ncN. (11) 

4) Capacity Constraints: To complete the characterization 
of the network, we point out that the data flows and the power 
allocations are coupled by the maximum aggregate rate that 
can be supported by the subcarriers of each link. In particular, 
the aggregate rate ^ik must not exceed the capacity of 

the k-th subcarrier of link £. 

To characterize the capacity constraints, we note that the 
transmission on link I G L and subcarrier fc G /C is 
composed of two parts. The first part accounts for the fraction 
of time over which this transmission is interference-free, 
whereas the second part accounts for the fraction of time over 
which this transmission interferes with other transmissions. 
To characterize the second part, we identify the interfer¬ 
ing links and the fraction of time over which these links 
are interfering. To do so, we note that, if subcarrier k is 
time-shared by links then the transmissions on 

links interfere with the transmission on link £ 1 . 

Hence, the SNIR expression for the transmission on link £i 
is -, where £'- denotes the index of the link 

connecting the node at which link £i originates to the node 
at which link £\ ends. Since links £\,...,£m are simultane- 
ously active on subcarrier k for a fraction of 7 )^ ^ , the 
expression for the data rate that can be communicated over 

link £1 is log2 (^1 -f ' Summing over¬ 

all possible combinations of the interfering links, the capacity 
constraint on the aggregate flow of link £i on subcatTier k can 
be expressed as 

E 4fk ^ log2(l + q£ik9Gk)+ 
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E E ^ 

m= 2 ^ 2 ...^mG£ 


,(fc) 


l0g2 1 + 




-i I 

1 + l^i=2 Itikgi'k 


( 12 ) 


B. Problem Formulation 

To ensure the feasibility of the rates generated by our design, 
the constraints in (l2li- (fT2T i must be satisfied. Combining these 
constraints yields the following optimization problem; 

max 

itvuJKid} 

subject to 

Routing constraints in (l2]i-(|4|i, 

Scheduling constraints in (|5]l-(|8]l, 

Power allocation constraints in (fTOl i and (fTTT i. 
Capacity constraints in (fT2T i. (13) 

The optimization problem in (fTTb is nonconvex because of the 
power allocation constraint in (fTTT i and the capacity constraints 
in (fT2l i. Examining (fOl l reveals that this problem shares some 
features with the GP standard form, cf. Appendix lA-ll To 
exploit this observation, in the next section we will perform a 
change of variables that will enable us to express the objective 
and all, but one set, of the constraints in a GP-compatible 
form. The residual constraints that do not comply with the 
GP standard form are approximated using the monomial 
approximation technique in Appendix IA-21 Under relatively 
mild conditions ||8l, iterative application of this technique is 
known to yield a solution of the KKT system corresponding 
to (fT^ . see e.g., fS), ifT^ . 

C. Generalized GP-Based Algorithm 

To cast (fTSl l in a form that is amenable to monomial 
approximation, we define two sets of variables, and 

which are related to and by the following 

maps: 

sW=log2fW, xfj =Wlog,r<fJ, 

n eAf\{d}, d evj e C, k £ 1C. (14) 

These maps are bijective, which renders recovering 

straightforward. The objective and the routing 
constraints in (fOl l can be cast in a GP-compatible 
form. In particular, the objective can be expressed as 

ridec n„eA^\{d} " and the routing constraints can 

be expressed as 

n n£^^\{dhd£v, ds) 

tGC keK 

>1, ieC,kelC,d£V, (16) 

> 1, n£Af\{d},deV. (17) 

The non-negativity constraints in Q and (fTOl i are inherently 
satisfied in the GP framework. The constraints in (|6]l and (fTTT i 
are already in a GP-compatible form. We now consider the 
constraints in (Q and ([8]l. The right hand side (RHS) of these 
constraints are zero, which makes them incompatible with 


the GP framework in Appendix lA-ll This problem can be 
alleviated by constraining their LHS to be less than an arbitrary 
small number e > 0, i.e.. 


^nl2 


(4S. + E E 


ii £ C,£2 £C\{h},k£ K, (18) 

/ 


+ E E 


m —3 £ 3 ---£rn^tC 


< e, 


£i e jC,£ 2 £C\{£i},k£ K. (19) 


The remaining constraints that are not GP-compatible are 
those in (fT2] |. Invoking the change of variables in (fTTl i. for 
£i £ C and k £ 1C, those constraints can be expressed as 


n ^etk < (1 + 


n n 


qiikgtik 

, . vm 

1 T /^i=2 Qtikgi'^k 



( 20 ) 


The RHS of (l20l i is amenable to the monomial approximation 
technique described in Appendix IA-21 171 . One approach to 
use this technique is to approximate all the terms in the 
RHS of (lini) by one monomial. This approach is overly 
complicated, and an alternative is to approximate each term 
by a monomial. The product of these monomials constitutes 
a monomial approximation of the RHS of (l20l . Hence, the 
constraint in (l20l can be approximated with 


n ^thk < x((i+) X 

deT) 

I / (fc) \ 

n n (21) 


where the functional Ad(■) is described in Appendix lA-21 Note 
that, { 7 ]^..^ } are variables and hence inseparable from the 
argument of M{-). 

Now, the problem in (fT3l can be approximated by the 
following GP: 


max 


n n ( 4 “’) 

dSi'D nGAf\{d} 



subject to 

Routing constraints in (fT5b-(fT7b. 

Scheduling constraints in (|6]l, (fTSl and ( fT^ . 
Power allocation constraints in (fTTl . 
Approximate capacity constraints in (fTTT i. 


( 22 ) 


Note that the relaxations in (fTSl and (O may result 
in infeasible subcarrier time-shares that do not satisfy the 
constraints in (|7|i and (O. To construct feasible schedules, the 
elements of T that are less than or equal to e are set to zero. 
Using a standard exponential transformation, the GP in (f22fl 
can be readily transformed into a convex optimization problem 
which can be solved in polynomial time using interior-point 
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TABLE II: Successive GP-based Algorithm for Solving 

1- Let = 0. Set accuracy to <5 > 0. 

2- Choose I and a feasible 

3- Solve the GP in I13> . Denote the value of the objective by U. 

4- While U - > <5, 

1/(0) ^ U, 

Solve the GP in 1131 . Denote the value of the objective by U, 
End. 

5- Remove the elements in F that ai‘e less than e. 

6- Use to recover and 


methods (IPMs) Q. This implies that (l22l i enables us to 
efficiently solve (fOl l approximately in the neighbourhood of 
any initial set 

Finding the global solution for the nonconvex problem 
in (fT3l) is difficult, whereas solving the approximated problem 
in (l22li is straightforward. To exploit this fact, we incorporate 
the formulation in (|2^ in an iterative algorithm, whereby the 
output of solving (l22l i for an initial point r(o)) 

is 

used as a starting point for the subsequent iteration. This 
technique is usually referred to as the single condensation 
method, e.g., mi, EO), and under relatively mild conditions, 
its convergence to a solution of the KKT system corresponding 
to (fOl) is guaranteed 18]. Since the original design problem 
is not convex, this system has multiple local solutions and 
the one to which the single condensation method converges 
depends on the initial point; some of the local solutions may 
be global ones. A summary of this algorithm is described in 
Table HI] 

In the next section, we will discuss a special case of this 
algorithm when time-sharing of subcarriers is not allowed. 
Before we do that, we now provide a brief discussion on 
the implementation of this algorithm. To begin with, we note 
that the algorithm in Table |II| is centralized, in the sense that 
the design is performed by a central entity that is aware of 
the network parameters. The signalling exchange between the 
nodes and the central entity, required to establish communi¬ 
cation in the considered framework, are described as follows. 
At the beginning of each signalling interval, the central entity 
prompts the nodes in the network to sequentially broadcast 
pilot signals of prescribed power levels. Subsequently, each 
node computes the subcarrier channel gains from all other 
nodes in the network. There is total of LK such gains, where 
L is the number of links and K is the number of subcarriers. 
Each node sends these gains along with its destination nodes, 
if any, and its priority weights to the central entity. The central 
entity performs the joint optimization of the power allocations, 
scheduling parameters and data routes as described in Table HIl 
It then forwards these decisions to all the nodes, possibly 
over a dedicated control channel. In particular, the information 
forwarded by the central entity include 1) the subcarrier index 
and the time allocated to each transmission. This information 
is provided by the set F. The cardinality of this set depends 
on the number of simultaneous transmissions allowed in each 
subcarrier, I. For instance, for 1 = 2, |r| = LK{L+ l)/2; 2) 
The power allocated to each transmission. This information is 


provided by the set {qik] and the cardinality of this set is LK-, 
and 3) The data rates at each transmitting and receiving node 
in the route of the stream intended for each destination. This 
information is provided in the set and the cardinality 

of this set is LKD, where D is the number of intended 
destinations. 


V. Special Case: Routing and Resource 
Allocation without Time-sharing 

In this section, we consider a design problem similar to 
the one described in Section IIVI but for the case when time¬ 
sharing of subcarriers is not allowed. This corresponds to the 
a special case in which the entries of F in Section IIV-A2I are 
restricted to be binary. This restriction results in a mixed inte¬ 
ger program which is generally difficult to solve. To overcome 
this difficulty, we capture the effect of the scheduling variables 
in the power allocation constraints. We will show that this 
approach will enable us to develop a design algorithm with 
polynomial-complexity. 

A. System Constraints 

1) Routing Constraints: These constraints are identical to 
those described in (|2]l-(|4|l. 

2) Power Allocation Constraints: In characterizing these 
constraints, we will use the method described in Section rV-A2l 
to denote the power allocated for transmission on subcarrier k 
of link £ by the variables {qek}- These variables must satisfy 
the non-negativity constraints in (fTOl i and the power budget 
constraint. These constraints, using ®, can be cast as 

EE atiqtk<Pni ncM. (23) 

keKeec 

Similar to the case considered in Section IIVI the nodes 
cannot simultaneously broadcast to multiple destinations on 
the same subcarrier. However, this requirement in the current 
case can be implicitly captured by the allocation of the link 
powers. In particular, for any subcarrier k € 1C, any node 
n € Af and any two links £ 1,^2 G C+{n), at least qi^k = 0 or 
qi^k = 0, i.e., 

ii, £2 € C, k € 1C, n € JV. (24) 

This constraint is significantly less involved than the one in (jTj. 

Similarly, the half-duplex requirement can be captured by 
ensuring that, for each node n G Af, if the power on subcarrier 
k of C+{n) is strictly positive, then the power allocated to 
this subcarrier on all the links in C-{n) is zero, and vice 
versa. Hence, the half-duplex requirement can enforced by the 
following constraints: 

^nei^ni 2 ^Ckqe 2 k = 0, £i, £2 c C, k € 1C, neAf. (25) 


Note that these constraints are simpler than their counterparts 
in (El). Also, note that (i24l) and (1251) are trivially satisfied if 
either link £i or £2 are not connected to node n. 

3) Capacity Constraints: In this case, the constraints in (fT^ 
can be readily seen to reduce to 


Ad) 

E d-l-^k 

w 

d 


< 


l0g2( 


1 + 


qiikgiik 


1 + J2t2GC\{ei} 9Ckge'2k 


(26) 










VI. Complexity Analysis 


B. Problem Formulation 

Using the characterization described in Section IV-AI the 
design problem can be cast as: 

max 'Y' 'Y' 

subject to 

Routing constraints in ©-dill, 

Power allocation constraints in (fTOl i. (l2^ - (l25l l. 

Capacity constraints in (l26l l. (27) 

The optimization problem in (l27l) is nonconvex because 
the RHS of (l26l) is the logarithm of a rational function, and 
therefore not concave. The equality constraints in (l24l) and (|25] | 
are not affine and hence, nonconvex. In the next section, 
we will develop a GP-based algorithm, analogous to the one 
described in Section ITV-CI to obtain a locally optimal solution. 


C. Proposed GP-based Algorithm 

The optimization problem in (l27l i. although nonconvex, 
is amenable to the GP-based monomial approximation in 
Appendix IA-21 To use this approximation, we use (fT4li to 
transform and {xg^} to {t^^} and respectively. 

Using these new variables, the routing constraints are readily 
expressed in a GP-compatible form as described in (fT5] l- (fT7] |. 

Substituting from (fT4li into (l26T l yields the following set of 
equivalent constraints: 

(l + ^ qi^kge'^k) n - 

foGAUi} 

fc e /C, £i G £, (28) 


which are, unfortunately not GP-compatible. Using the mono¬ 
mial approximation technique in Appendix IA-21 yields the 
following approximation of (l28T l in the neighbourhood of 


(i+^ n ’’Gfc - n {^Gk/qfjk) 

Ger;\Ui} d&v faez; 


kelCJi^C, (29) 


where is the initial power allocation, = 1 + 

qefkde'.k, and = qf^St'-^k!^Gk- 
Analogous to the case considered in Section IIV-CI (l24l i 
and ( |25] | are replaced with the GP-compatible inequality con¬ 
straints. The joint design of data routes and power allocations 
in (l27l i can be approximated with the following GP: 


max 


n n 

dGU n£AZ\{d} 


„(<!) 


subject to 

Routing constraints in (fTSTl-lfTTIi. 

Power allocation constraints in (l2^ - (l25l l (relaxed versions), 
Approximate capacity constraints in ( |29] |. (30) 


A locally optimal solution of ( |27] | can be obtained by solv¬ 
ing ( l30b iteratively using the single condensation method 
described in Section IIV-CI 


In this section we examine the computational complexity 
required for solving the problems described in Sections IIV-CI 
and IV-CI for the cases with and without time-sharing, re¬ 
spectively. The algorithms in these sections iteratively solve 
the families of the optimization problems in (l22l i and (l30l) . 
Being in a GP-compatible form, these problems can be readily 
converted into convex forms and can be efficiently solved 
using IPM-based solvers. 

In IPM, the objective and inequality constraints are used 
to construct a log-barrier function which is minimized along 
a central path using Newton’s method. The complexity of 
each Newton step grows with the cube of the number of 
inequality constraints and the number of Newton steps can 
be bounded if the log-barrier function is self-concordant ll2Tll . 
cf. Appendix |B] In that case, the number of Newton steps 
can be shown to grow with the square root of the number of 
inequality constraints 11211 . 

Unfortunately, the log-barrier functions related to the prob¬ 
lems in (I 22 I 1 and (|30] | are not self-concordant. To circumvent 
this difficulty, we introduce a set of auxiliary variables and 
constraints which, although redundant, enables us to construct 
self-concordant log-barrier functions. Using these functions 
and the results in ||2T], we arrive at the following proposition: 

Proposition 1. The complexity of solving (l22l) with IPM-based 
solvers is of order 

O ^ ^LKN + N + DiN - 1) + 2KY2 

and the complexity of solving dlOll with IPM-based solvers is 
of order 

O (^{LK{3L + 2) + N + D{N - 1))^'®) . 

Proof: See Appendix iBl ■ 

The first statement of Proposition [T] pertains to the general 
case with time-sharing and frequency-reuse. This statement 
shows that the complexity of solving the problem in (l22l) 
is polynomial in L for small values of I. The complexity 
of solving (|2^ can be further reduced by combining the 
broadcasting constraint and the half-duplex constraint in (|7]i 
and ®, respectively. In particular, examining these constraints 
reveals that they are related to the network topology and do not 
depend on the channel conditions. Hence, these two constraints 
can be enforced by pruning the set T prior to solving (fT^ or its 
approximated version in (l22l i. The pruning rule is as follows: 
For each I and P G C, if either = 0 or = 0, 

the corresponding time-shares in © and © are removed from 
the set r. Unfortunately, we have not been able to obtain a 
closed form of the cardinality of the resulting T. However, the 
reduction in complexity, at least for small networks, appears 
to be significant. For instance, for fully connected networks 
with = 4 nodes and L = N{N — 1) = 12 links, |r| is 
reduced from 4095 to 40. 

The second statement of Proposition[T]pertains to the special 
case in which time-sharing is not allowed. This statement 
shows that the complexity of solving (l30t is polynomial in the 
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TABLE III: Normalized Channel Gains, {gik\, in Example 1 [dB], 



link 1 

link 2 

link 3 

link 4 

link 5 

link 6 

subcanier 1 

-6.1 

-11.1 

1.86 

-12.3 

-13.5 

42.4 

subcarrier 2 

-5.7 

1.57 

5.64 

-5.98 

-19.2 

40.7 


link 7 

link 8 

link 9 

link 10 

link 11 

link 12 

subcanier 1 

15.82 

-2.43 

-5.82 

-6.9 

35.0 

3.304 

subcanier 2 

-3.38 

-5.61 

0.871 

-0.6 

41.2 

-14.5 


number of nodes, N, and the number of subcarriers, K. In par¬ 
ticular, it grows as . Another case in which the design 

complexity is polynomial is the one in which the subcarriers 
are time-shared but not frequency-reused 0. In that case the 
design complexity 0{{LK{A+D)+N-\-K-\-D[N — 

Hence the special cases with either no frequency-reuse or no 
time-sharing have polynomial complexity. 

VII. Simulation Results 

In this section we provide numerical results to evaluate 
the performance of joint routing and resource allocation al¬ 
gorithms for the cases with and without time-sharing. 

The locations of the nodes are randomly generated and 
evenly distributed over a 300 x 300 m^ square. The nodes are 
assumed to have identical power budgets, i.e, = P, Vn G 
M, and the available frequency-selective channel is partitioned 
into a set of frequency-flat Rayleigh fading channels with 
the values of pathloss (PL) and shadowing components ob¬ 
tained from the non line-of-sight communication of indoor 
hotspot (InH) scenario in the IMT-Advanced document 
According to ll22ll . the PL component on link £ G £ is given 
by 

PL = 43.31ogio(df) + 11.5 + 201ogio(/c), (31) 

where di is the length of link I in meters and /c is the 
carrier frequency in Gigahertz which, in our simulations, is 
set to /c = 3.4 GHz. The shadowing component is assumed 
to be log-normal distributed with a mean of 0 dB and a 
standard deviation of 4 dB. The Rayleigh fading component 
is generated by the envelope of a zero-mean unit-variance 
complex Gaussian-distributed random variable. The available 
bandwidth around each subcarrier is set to W = 200 KHz 
and the noise power density at receivers is set to Nq = 
-174 dBm/Hz. 

The results reported herein are obtained using the CVX 
package ll23l with an underlying MOSEK solver ll24ll . The 
value of e in (|7]l and ® is set to 10“^. 

Example 1: (Joint Routing and Resource Allocation with 
Timesharing) Consider an exemplary network with N = A 
nodes. In this network, nodes 3 and 4 wish to communicate 
with nodes 2 and 1, respectively, over K = 2 subcarriers. In 
particular, for destination node d = 1, the source is node n = 4 
and nodes {2, 3} are potential relays, and, for destination node 
d = 2, the source is node n = 3 and nodes {1,4} are potential 
relays. The considered network has L = 12 directional links 
and therefore the channel matrix has 12 x 2 elements. The 
channels are assumed to be static and their normalized gain in 
dB, i.e., 10 log^o is given in Table Hill In this example, the 
power budget of each node is set to P = 20 dBm, the number 
of simultaneous transmissions is set to / = 3 and the two 


TABLE IV: Power Allocations (mW) in Example 1. 


n = 1 

n = 2 

n = 3 

n = 4 

93,2 = 23 

94,2 = 56.5 

97.1 = 6.5 

99.2 = 13.5 

910,2 = 24.5 

911.1 = 0-3 

912.1 = 16.5 


rates, and 54^^ are assigned equal weights, i.e., = 

= 1. Since in this example time-sharing is allowed, the 
algorithm in Section IIV-CI is used to generate the data routes, 
time-sharing schedules and power allocations. 

The sum-rate yielded by the algorithm in Section IIV-CI is 
7.4 b/s/Hz. The data routes generated by this algorithm are 
illustrated in Figure [3 For ease of exposition, the network 
in this example is split into the two sub-networks: the one 
in Figure |3(a)| depicts the routes of the data intended for 
destination d = 1, and the one in Figure |3(b)| depicts the 
routes of the data intended for destination d = 2. The complete 
network is the superposition of the two sub-networks. For 
instance, the data transmitted over link 7, connecting node 
3 to node 1, is 4.8 b/s/Hz, of which 2 b/s/Hz is intended for 
destination d = 1 and 2.8 b/s/Hz is intended for destination 
d = 2. 

The time-sharing schedules of the subcarriers generated by 
the algorithm in Section IIV-CI are provided in Figure |4] It can 
be seen from this figure that subcarrier k = I is both reused 
and time-shared, whereas subcarrier fc = 2 is only time-shared. 
Figures [3 and |4] imply that link 7, connecting node 3 to node 
1, and link 11, connecting node 4 to node 2, carry the data 
intended for both destinations on the same subcaiTier, k = 1, 
during the same time interval. The fact that our designs enforce 
half-duplex requirement can be inferred from these figures. 
For instance. Figure |3(a)| shows that node 3 uses the same 
subcarrier, k = 1, for its transmission and reception on links 
7 and 12, respectively, but Figure |4] shows that transmission 
and reception occur during different time intervals. 

The power allocations yielded by the algorithm in Sec¬ 
tion ITV^ are shown in Table El This table shows that, 
because of frequency-reuse, the nodes do not necessarily use 
their total power budgets. This is due to the fact that, in 
this scenario, when a node increases its transmission power, 
it inflicts high interference on other transmissions. This is 
in contrast with the situation considered in 0, wherein 
frequency-reuse is not allowed and increasing the transmitted 
power of a node does not affect the transmissions of the other 
nodes in the network. □ 

Example 2: (Joint Routing and Resource Allocation without 
Timesharing) Consider an exemplary network with N = 6 
nodes. In this network, as before, nodes 3 and 4 wish to 
communicate with nodes 2 and 1, respectively, over AT = 4 
subcarriers. In particular, for destination node <7=1, the 
source is node n = 4 and the other nodes, i.e., {2,3,4, 5} are 
potential relays, and, for destination node d = 2, the source is 
node n = 3 and nodes {1,4, 5, 6} are potential relays. 

The considered network has L = 30 links and therefore the 
channel matrix has 30 x 4 elements. For space considerations, 
this matrix is not provided, but since the channel gain on each 
subcarrier is dominated by the PL component, we provide 
the coordinates of the nodes in the 300 x 300 m^ square; 
calculating the PL components from these coordinates is 
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\ 4.8 b/s/Hz 


^ode^ 


^^ode^ 


= 4.8 b/s/Hz 


(a) 



Fig. 3; Data routes forint/ = 1, |(b)|d = 2 in Example 1. 


k = 1 


Link 11 
Link 12 


k = 2 


Link 3 
Link 4 
Link 9 
Link 10 


Fig. 4: Time-sharing schedules of the subcarriers in Example 1. 

straightforward, cf. OTI ). The coordinates of the nodes are 
{(283,202), (191,208), (287,20), (72,76), (201,67), (86,200)} 

Setting the node power budgets to P = 20 dBm and 
assuming that both rates have equal weights, = 1, 

the joint design algorithm in Section IV-CI yields a sum-rate of 
9.1 b/s/Hz. The data routes and power allocations obtained by 
this algorithm are shown in Figure|5]and Table [Vl respectively. 
For instance, in Figure |5] subcarrier fc = 1 is shown to be used 
twice and due to the half-duplex constraint, transmission and 
reception take place over distinct subcarriers at each node. We 
will later show the advantage of the proposed algorithm over 
the algorithms in which frequency-reuse is not considered. □ 

Example 3: (Average Weighted-Sum Rate Comparison) In 
this example, we use Monte Carlo simulations to evaluate the 
average performance of the joint designs with and without 
time-sharing when the channels are time-varying rather than 


TABLE V: Power Allocations (mW) in Example 2. 


n = 1 

n = 3 

72 = 4 

72 = 5 

91,1 = 12 

911,2 = 45 
914,1 = 55 

916,4 = 100 

922,3 = 25 


Fig. 5; Data routes for |(a )| d = l,|(b)|d = 2 in Example 2. 


Static as in Examples 1 and 2. We consider a network with 
N = 4: nodes in which nodes 3 and 4 wish to communicate 
with nodes 2 and 1, respectively, over K = A subcarriers. The 
number of simultaneous transmissions is set to / = 3 and the 
simulation results are averaged over 10 independent network 
realizations. 

The average weighted-sum rates yielded by the algorithms 

in Sections IV-CI and IIV-CI for the values of P ranging from 0 

to 30 dBm are depicted in Figures |6(a)| and |6(b)| for the cases 
( 2 ) ( 1 ) ( 2 ) ( 1 ) ’ ’ 
of w\ = and Wg = w\ , respectively. These figures 

also provide a comparison with the weighted-sum rates yielded 

by the designs in which frequency-reuse is not considered 0. 

As can be seen from Figure|6] the weighted-sum rate yielded 
by the joint design with both time-sharing and frequency- 
reuse outperforms the designs in which either time-sharing or 
frequency-reuse is exclusively considered, but at the expense 
of increased complexity. For instance. Figure |6(b)| suggests 
that, at the sum-rate of 12 b/s/Hz, the proposed design with 
both time-sharing and frequency-reuse yields a power advan¬ 
tage of 4 dBm over the designs in which either time-sharing 
or frequency-reuse is exclusively considered and a power 
advantage of 8 dBm over the design in which neither of these 
techniques is considered. This figure also suggests that, for 
values of P less than 15 dBm, the design with frequency-reuse 
but without time-sharing yields better performance than the 
design with time-sharing but without frequency-reuse in 0. 
However, for values of P higher than 15 dBm, the design with 
time-sharing but without frequency-reuse performs better than 
the one with frequency-reuse but without time-sharing. This 
phenomenon can be attributed to the effect of interference. At 
low powers, the effect of interference is small and frequency- 
reuse performs generally better than time-sharing. In contrast, 
at high powers, the effect of interference is more severe 
and time-sharing performs generally better than frequency- 
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Fig. 6; Average weighted-sum rate comparison for |(a)| w^^ = 
and |(b)| w^^^ = 

reuse. As expected, the design with neither time-sharing nor 
frequency-reuse has inferior performance. □ 

Example 4: (Joint Routing and Resource Allocation: Gen¬ 
eralized Algorithm) In this example, we evaluate the per¬ 
formance of the algorithm developed in Section IIV-CI We 
consider a snapshot of a network with N — 5 nodes and 
L = 20 links (links with a distance more than 150 m are 
neglected). In this network nodes 3 and 4 wish to communicate 
with nodes 2 and 1, respectively, over K = 4 subcarriers. The 
number of simultaneous transmissions is set to / = 20, / = 3 
and 1 = 2, which results in T with 2^° — 1, 190 and 110 
variables, respectively. 

The sum-rate yielded by the generalized algorithm with 
different values of I is depicted in Figure |2l For comparison, 
this hgure also shows the rates yielded by the special case 
in Section [V] As can be seen from Figure |7] the algorithm 
with 1 = 2 and 3 yields rates that are slightly less than 
the rate yielded by the algorithm with I = L, however 
with a signihcantly less computational complexity. In fact, 
the complexity of the algorithm with I L is polynomial, 
whereas that of the algorithm with I = L is exponential in L. 
This feature renders the algorithm with / < 3 more attractive 
for designing large networks with potentially rapid channel 
variations. From Figure [T] it can be seen that the gap between 
the rates yielded by different values of I decreases as the 
power budget increases. This is because as power increases, 
interference becomes more severe, which causes the reuse 
of a particular subcarrier on multiple links less benehcial. 
It can be also seen from this hgure that, most of the gain 
of frequency-reuse is mustered by only considering two or 


three simultaneous transmissions, i.e., / < 3. This implies that 
increasing I trades complexity for performance. In particular, 
as I increases, the performance of the algorithm becomes 
closer to that of the one with I = L, but at the expense of 
increased comnlexitv. 



Fig. 7: Sum-rate generated by tbe generalized algorithm for different 
values of I. 

In Figure |8(a)| we investigate the convergence behaviour 
of the generalized algorithm. We consider an instance of a 
network in which the power budget of all nodes is set to 
P = 10 dBm. It can be seen from this hgure that, in addition 
to being signihcantly less computationally demanding, the 
algorithm with a lower value of I exhibits considerably faster 
convergence than that of the one with a higher value of I. This 
convergence can be further ameliorated by choosing the initial 
point more carefully, for instance, by choosing this point to 
be the solution yielded by algorithm in Q for the case with 
no frequency-reuse. 

To illustrate the effect of random initialization of the gen¬ 
eralized algorithm, in Figure |8(b)| the value of the objective 
to which the generalized algorithm with 1 = 2 converged 
is shown for 80 random instances of feasible initial points, 
£ [0,P]^^ X [0,1]I^L It can be seen from this 
hgure that although the algorithm is relatively sensitive to the 
choice of the initial point, hnding initial points that result in 
‘good’ local maxima is generally easy. □ 

Example 5: (Average Rate-Region Comparison) In this 
example we provide the rate regions that can be achieved by 
the algorithms in Sections I V-CI and llV-CI when P = 10 dBm. 
These regions are obtained by varying the weights (s4^\ 
over the unit simplex, i.e., > 0, > 

0, ^4^^ = l}. are depicted in Figure|9] A compar¬ 

ison between these rate regions and the ones corresponding 
to the case when frequency-reuse is not considered El is 
also provided in this hgure. As can be seen from Figure |9] 
the rate region corresponding to the design with both time¬ 
sharing and frequency-reuse properly contains the rate regions 
corresponding to the designs in which either time-sharing or 
frequency-reuse is exclusively used. It can be also seen that 
restricting the number of simultaneous transmissions to be less 
than three suffices to achieve most of the frequency-reuse gain 
and with less computational complexity. □ 

VIII. Conclusion 

In this paper we focused on the joint optimization of 
data routes, subcarrier schedules and power allocation in a 
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i 7.4 

tn 
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Number of experiments 

(b) 

Fig. 8: |(a)| Convergence behaviour and |(b)| performance of the 
generalized algorithm with different initial points. 




half-duplex multicarrier network when each subcarrier can 
be reused by multiple links. The goal is to maximize a 
weighted-sum of the rates communicated over the network. 
The considered network is generic in the sense that it sub¬ 
sumes many structures including cellular and device-to-device 
communications as special cases. We considered two instances 
of this problem: 1) when each subcarrier can be time-shared 
by multiple links; and 2) when time-sharing is not allowed and 
a subcatrier, once assigned to a set of links, will be used by 
those links throughout the signalling interval. The joint design 
in the first instance results in superior performance but with 
high complexity. The second instance is a special case of the 
first one and can be parameterized using a significantly smaller 
number of variables. 

The joint design problem in both instances is nonconvex 
and locally optimal solutions are obtained using a GP-based 
monomial approximation technique. Numerical results show 
that the designs developed in both instances yield performance 
that is significantly better than that of their counterparts in 
which frequency-reuse is not allowed. 


Appendix A 

The GP Standard Form and Monomial 
Approximation 

1) The GP Standard Form: For self-containment, in this 
appendix we will review the standard GP form. A GP opti¬ 
mization problem can be readily transformed to an efficiently 
solvable convex one. To provide the standard form of a GP, 
let z € M” be a vector of positive entries. A monomial 
in z is defined to be a function of the form cq Hi 

a posynomial in z is defined to be a function of the form 
E/=r Cj nr=i where Cj > 0, {a,} and {a^}, are 
arbitrary constants, j = 0,1,... , J, and i = l,...,n. A 
standard GP I?], lITTl is an optimization of the form: 

min /o(z), 

2 

subject to fi{z) <1, i = 1,... ,m, (32) 

5i(z) = 1, i = l,...,p, 

where {fi} are posynomials and {gi} are monomials. 

2) Monomial Approximation: A monomial approximation 
of a differentiable function h{z) > 0 near is given by 

its first order Taylor expansion in the logarithmic domain ||6l, 

( 0 ) 

lU. Defining Pi = 7ipoyy|i7L=^(o), we have M[h{z)) = 
( \ 

,where Ad(-) is the monomial approxi¬ 
mation. This approximation will be used to provide local GP 
approximations in the neighbourhood of a given initial point. 


Appendix B 

Proof of Proposition[T] 

For simplicity, we will begin by proving the second state¬ 
ment of Proposition [T] In this proof, we will show that, by 
including redundant constraints, the log-barrier function of the 
problems in (l30l l and (l22li can be cast in a self-concordant 
form, which has the following definition lITTI : 


Definition 1. A function f 
concordant if, for all a:, u € 1 


s e 


R is said to be self- 
such that X sv is in 


the domain of f and 


§^f{x-\-sv) < 2^/(a;-f □ 


A. Proof of the Second Statement of Proposition Q] 

To determine the complexity of solving the problem in (l30l) . 
we begin by converting this problem into a convex one. Using 
standard exponential transformations, we write 

tn'^ = exp , n&M\ {d}, d G D, 

(d) 

rff} = exp(^ln(2)^-^, £ G C, k G IC, d € V, 

yik = exp(g^fc), £ G C, k G fC. (33) 

Substituting the variables in (l30l l with the ones in (l33]) and 
taking the logarithm of the obtained objective and constraints 
result in a convex optimization which can be solved effi¬ 
ciently using the IPM technique. To use this technique, a 
log-barrier function is synthesized from the objective and 
inequality constraints. The complexity analysis of the IPM 
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technique is simplified when the log-barrier function is self- 
concordant im, cf., Definition [T] The log-barrier function 
corresponding to the convex form of ( |30] | can be written as 

(j) =+-1^, (34) 

n d 

where ip represents the component of the log-barrier function 
associated with the inequality constraints in the convex form 
of (l30t . To examine whether (p is self-concordant, we note that 
the converted objective and the inequality constraints corre¬ 
sponding to ( fTfil l. ([TtI i and the relaxed versions of (l24l i and (l25T l 
are linear and therefore their corresponding components in the 
log-barrier function are self-concordant ED. Hence it remains 
to consider the self-concordance for the constraints in (l29l l 
and (|2^ . For simplicity, we write the posynomial constraint 
in (|2^ in the standard form in (l32l i. After changing the 
variables and taking the logarithm of both sides, this constraint 
can be written in a general form as 

log(^ exp(aiai -f bi^i + a)) < 0, (35) 

i 

where {ai],{l3i} are the optimization variables and 
{tti}, {bi}, {ci} are constants. The component corresponding 
to the constraint in ( iTSl l in the log-barrier function can now 
be expressed as 

- log(- log ^ exp(aiai -f bi/3i -f Ci)). (36) 

i 

To ensure that ( l36l l is self-concordant, we introduce auxiliary 
variables, Xi, to bound the exponentially transformed variables 
in (ITST i. Using these new variables, the constraint in dTSl l can 
be replaced with the following set of constraints m- 

< 1 , 

i 

x^>o, 

UiUi + bif3i +Ci- log Xi < 0. (37) 

Now the associated log-barrier function of the constraints 
in (I37li can be shown to be self-concordant, cf., lEI] Ex¬ 
ample 9.8]. For the constraints in (l23t . we follow the steps 
analogous to the ones used with the constraints in ( |29] |. 
In particular, by introducing new auxiliary variables, we 
construct a self-concordant log-barrier function. Using this 
function, the complexity can be shown to be proportional 
to where m is the number of inequality constraints. 

Hence, the complexity of solving (l30l l can be bounded by 
O (^{LK{3L + 2) + N + D{N — 1))^'^^, which completes 
the proof of the second statement of Proposition [T] 

B. Proof of the First Statement of Proposition Q] 

The proof of the first statement of Proposition [T] follows 
from arguments similar to the one used in the proof of the 
second statement and is omitted for brevity. For the first 
statement, the number of inequality constraints can be readily 
verified to be 2LKN + N + D{N -l) + 2KY!i=i {^), which 
yields the first statement of Proposition [T] 
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